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We show how fragmentation of a Bose-Einstein condensate can occur given repulsive inter-particle 
interactions and a non-uniform external potential. 

I. INTRODUCTION 

It is customary to approximate the iV-body ground state of a system of spinless bosons by the Hartree-Fock state 
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where 4>a{r) is a normalized single particle wavefunction jjj]. This describes an accumulation of N particles in the 
wavefunction (j)o(r). One can, however, also consider states of the form 
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where 4>i{r) and 4> 2 {r) are normalized and orthogonal single-particle wavefunctions, and where Ni + N 2 = N. Such a 
state describes an accumulation of Ni particles in the wavefunction 4>i(r), and N 2 particles in (j> 2 (r). We distinguish 
states of the form \N, </>o) from those of the form \Ni, </>i ; N 2l 4*2) by designating the first as 'single condensates' and the 
second as 'dual condensates'. It is evident that one can also introduce states describing the accumulation of particles 
in an arbitrary number of orthogonal single-particle wavefunctions. Such states can be referred to collectively as 
'fragmented condensates' In this paper, we show how a non-uniform trap potential can encourage fragmentation. 
In particular, we demonstrate that in a double minimum potential a dual condensate can have an expectation value 
of the energy smaller than that of a single condensate, and therefore can constitute a better approximation to the 
fully interacting ground state and a better starting point for treatments beyond Hartree-Fock. 

This result may be surprising to some and obvious to others. Specifically, it seems remarkable in light of Nozieres' 
argument against fragmentation of condensates [||. On the other hand, it seems natural to describe bosons in well- 
separated traps by a product of single condensate states for each trap; for two traps, such a product state would be a 
dual condensate wherein each of <^i(r) and 4> 2 (r) is centered in only one of the traps. Indeed, this has been assumed 
implicitly in the literature when discussing the interference of condensates ||. Our analysisjustifi.es this intuition and 
specifies the circumstances under which Nozieres' argument is inapplicable. 



II. DESIDERATA FOR FRAGMENTATION 



We consider a two-particle interaction that is local and completely repulsive: V(r, r') = g5(r — r'), The interaction 
strength is given by g — Aira sc Ti 2 /m, where a sc is the s-wave scattering length, assumed to be positive, and where m 
is the mass of the particle. For simplicity we restrict ourselves to wavefunctions <po(r), <fii(r), and 4> 2 (r) that are real. 
Consider the expectation value of the energy for a single condensate \N, </>o) , 
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where e(4>) — J (j)(r) — ^;V 2 + U(r) c/)(r)d 3 r is the single particle energy associated with </>(r), and U(r) is the trap 
potential. For the dual condensate jiVi^i; N 2 , cf> 2 ) , 
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If we assume that the particle densities are similar throughout the single and dual condensates, Ai0 2 (r) + iV202(r) ~ 
N(()q(y), then the pieces of the interaction energies that scale as N 2 are 

Er ad = \gN 2 J <j> 4 (Y)d 3 Y 
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We can think of this in the following way: while the Hartree contributions to the interaction energy of the single and 
dual condensates are essentially equal, the interaction energy of the dual condensate also includes an exchange term, 
proportional to N1N2, arising from symmetrization. 

This extra exchange energy was identified by Nozieres and is the basis of his argument against fragmentation 
of the condensate, since in the limit of a macroscopic number of particles, £^ uad exceeds £^ uad by a macroscopic 
energy. However, this argument assumes that the wavefunctions have a similar spatial extent (j>\(r) ~ 0§(r) ~ 0o( r )- 
If instead one considers a dual condensate wherein 0i and 2 have little density overlap, then the exchange term may 
be small compared to other terms in (^) that are linear in N, and we must consider the relative weights of all the 
terms. 

The energies £^ uad and E£ uad overestimate the interaction energy since they incorporate a particle's interaction 
with itself. Thus, to obtain the proper interaction energy, we must add to _E^ uad a term linear in N 

E^ = ~gN J 4(r)<Pr. (5) 

Similarly, we must add to £;^ uad the linear term 

= ~\gN t J 0?(r)d 3 r - l -gN 2 J 0|(r)<i 3 r. (6) 

If 0i and 2 are more localized than (f> , then / 0f 2 (r)d 3 r > J </$(r)d 3 r. It follows that E l j n -E£ n < ~\gN J o l (r)d 3 r. 
Thus, the linear piece of the interaction energy favors the dual condensate. 

Admittedly, it will cost more single particle energy to place particles in localized wavefunctions like 0i and 02 than 
in 0o, so the single particle energy favors the single condensate. Thus we see that there is a competition between the 
interaction energy benefit that may arise from separating the bosons, and the single particle energy cost of establishing 
this separation. The role of the trap potential in fragmentation is now clear: for a non-uniform potential, one can 
construct wavefunctions that are confined to the valleys of the potential and which have little density overlap without 
a large cost of single particle energy. A non-uniform potential therefore encourages spatial fragmentation. 

III. FRAGMENTATION IN A DOUBLE WELL POTENTIAL 

We now turn to a simple model of the condensate in a potential well that is symmetric along each Cartesian axis, 
but which along one axis (x) exhibits a double minimum. We wish to consider which of the single or dual condensate is 
energetically favored as a function of the strength of the central barrier of the potential. As is customary, one minimizes 
the energy of the single condensate with respect to variations in 0o to obtain a non- linear Schrodinger equation (NLSE) 
for 0o. Due to the symmetry of the potential, we expect that the energetically favored dual condensate will have the 
particles distributed equally between wavefunctions that are mirror images of one another across x = 0, that is, 
0i (— x,y,z) = cf>2(x,y,z) and N\ = N2 = N/2. Assuming a dual condensate of this form, we can minimize its energy 
with respect to variations in 0i and 02 to obtain two coupled NLSEs for 0i and 02- 

Let us first consider the case of an infinitely strong barrier; all wavefunctions will then have zero amplitude at the 
center of the barrier. The single condensate that minimizes the expectation value of the energy has a degeneracy 
associated with the different possible symmetries of the wavefunction across the barrier. Let us denote the symmetric 
solution of the single condensate NLSE as S and the antisymmetric solution as a . We can then construct two 
wavefunctions 0/ = 2 _1 / 2 (0 S + a ) and r = 2~ 1 / 2 (0 S — Q ) that are confined respectively to the left (x < 0) and 
right (x > 0) of the barrier, and which have the same single particle energy as S . 

In the single condensate \N, S ), the magnitude of the interaction energy is ^gN(N — 1) J0 2 (r)<i 3 r. This can be 
understood as follows: there are N particles interacting with N—l others with an interaction strength of g J 2 (r)d 3 r, 
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and the factor of 1/2 is due to over-counting. Now consider the particular dual condensate that is of the form 
| AT/2, fa, N/2, fa) , noting that this state is not necessarily the solution of the dual condensate variational problem. The 
interaction energy of this state originates solely from the interaction of particles with others in the same wavefunction, 
since the density overlap of the two wavefunctions, J (f>f(r)fa 2 (r)d 3 r, is zero. Particles in this dual condensate interact 
twice as intensely with one another than in the single condensate since J (j)f(r)d 3 r = J fa(r)d 3 r = 2 j fa(r)d 3 r. 
However, each of the ./V particles sees only N/2 — 1 others rather than N — 1 others. Including the over-counting 
factor of 1/2, the total interaction energy of this dual condensate is ^gN(N — 2) J fa(r)d 3 r. This interaction energy 
is clearly less than that of the single condensate. Since the single particle energies are the same for the two states, it 
follows that the dual condensate \N/2, fa, N/2, fa) is energetically favored overall. Further, the dual condensate that 
is obtained from solving the variational problem may have an energy that is lower still. 

For small barrier strengths we expect the single condensate to be energetically favored. This is because a dual 
condensate that has little density overlap of fa and fa will necessarily have a greater single particle energy than 
the single condensate. Alternatively, if the single particle energies arc comparable, the wavefunctions fa and fa will 
necessarily have a large density overlap. It is conceivable that if the interaction strength were strong enough, it could 
become energetically favorable for the condensate to fragment even in the absence of any trap potential. However, if 
the scattering length, a sc is much less than the average width of the ground state of the trap, it can be verified that 
a single condensate is preferred at small barrier strengths. 

To determine precisely the barrier strength at which the dual condensate becomes energetically preferred over the 
single condensate, one must solve the variational problems for both types of condensate and compare the energies 
obtained. 



IV. THE EXPERIMENTAL SIGNATURE OF FRAGMENTATION 



The most conspicuous difference between spatially fragmented and non-fragmented condensates is the extent of 
spatial coherence, or long range order, in the system. Specifically, let us imagine a double well potential with a large 
central barrier and ask how one could distinguish a single condensate | N, 4> s ) where <p s is a doubly peaked symmetric 
wavefunction, from a dual condensate | N/2, fa, N/2, fa) where 4>i and fa are mirror images of one another centered 
on the left and right sides of the barrier respectively. The answer to this question can be found in the recent literature 
P); essentially, the difference is revealed through interference experiments. 

A single run of such an experiment proceeds as follows: after the condensate is formed, the trap potential is removed 
instantaneously and the particles are allowed to expand freely until the wavefunctions from either side of the barrier 
have overlapped, at which time the positions of the particles are measured. Assuming the particles do not interact 
much during the expansion, the wavefunctions for each particle will evolve independently. If one models the particle 
detection in a manner analogous to the standard theory of photon detection, as is proposed in then both the 
single and dual condensates will generate an interference pattern. The single condensate exhibits interference since 
each particle interferes with itself just as it would in a double slit experiment. In the dual condensate, the probability 
density for each particle is essentially uniform before the first detection, but as is described in || , the first detection 
modifies the state such that the probability densities for the second and subsequent detections are no longer uniform, 
and an interference pattern emerges. 

The difference between the single and dual condensates can only be seen by a comparison of how the spatial phase 
of the interference pattern varies from one experimental run to the next. For the single condensate, the spatial phase 
is the same in every run, while for the dual condensate, it varies randomly. 

In the experiment of Andrews et al. Q a cloud of sodium atoms was cooled to below the phase-space density 
threshold for Bose-Einstein condensation in a double well potential. Upon releasing and detecting these atoms, an 
interference pattern was observed. Whether the spatial phase of interference patterns from different experimental runs 
are fixed relative to one another or not has not yet been determined due to mechanical instabilities in the apparatus. 
Although it was not our intention in this paper to provide a detailed model of this experiment, it is clear that the 
considerations outlined herein will be relevant to any such model. 



[1] Another common approach is to associate with the condensate a coherent state over Hartree-Fock states, \a,<j>o) = 
exp(a J d 3 r0o(r) , pt(r)) |vac) , where a is a complex parameter with modulus squared equal to N. It is unclear how one 
could describe the effects discussed herein using such an approach. 
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